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Based on a sample of 1.31 10° J /w events collected with the BESIII detector, the matrix elements 
for the decays n — nnr and n > nr°n® are determined using 351,016 n > (n > yy)at a7 
and 56,249 n! — (n > yy)r°r® events with background levels less than 1%. Two commonly used 
representations are used to describe the Dalitz plot density. We find that an assumption of a linear 
amplitude does not describe the data well. A small deviation of the obtained matrix elements 
between n — natn and n! > mgax?«? is probably caused by the mass difference between charged 
and neutral pions or radiative corrections. No cusp structure in n — rn^? is observed. 


PACS numbers: 13.66.Bc, 14.40.Be 


I. INTRODUCTION 


The 7 meson is well established and its main de- 
cay modes are fairly well known [1]. However, ņ de- 
cay dynamics remains a subject of extensive theoreti- 
cal studies aiming at extensions of the chiral perturba- 
tion theory (ChPT). The two dominant hadronic decays, 
n! + nxt and n” — nr°r® (called charged decay 
mode and neutral decay mode throughout the text, re- 
spectively), are believed to be an ideal place to study a7 
and nm scattering [2], which may lead to a variation in 
the density of the Dalitz plot. Several extensions of the 
ChPT framework have been applied to investigate the 
Dalitz plot of 7’ — nx [3-6]. 

In experimental analyses, the Dalitz plot for the 
charged decay mode is usually described by the following 
two variables: 


2m; T, 
AURE TERE y= nt ee. (1) 

Q m. è Q 
For the neutral decay mode, the Dalitz plot has twofold 
symmetry due to the two 7°s in the final state. Hence, 
the variable X is replaced by 


x= V3|Tzo =Tal 
Q 


Here, Tp and T, denote the kinetic energies of pion and 7 
in the 7/ rest frame, Q = m, — my — 2Mr, and Mz, Mn, 
and m, are the masses of pion, 7, and 7’, respectively. 
Generally, the decay amplitude squared is parameterized 
as 


(2) 


|M(X, Y)? 2 N(1-- aY -- bY? +cX +dX7+...), 
(3) 
which is the so-called general representation. Here a, b, 
c and d are free parameters and N is a normalization 
factor. The terms with odd powers in X are forbidden 
due to charge conjugation symmetry in n — nrn and 
the wave function symmetry in n’ > n7°7r°. By consid- 
ering the isospin symmetry, the matrix elements for the 
charged and neutral decay modes should be the same. 
However, a small discrepancy, observed in previous mea- 
surements [1], is expected due to the mass difference be- 
tween the charged and neutral pion, or due to radiative 
corrections for the n + nnn mode [7]. 
A second parameterization for the decay amplitude 
squared used by previous experiments assumes a linear 


amplitude in Y and keeps the polynomial expansion in 
X 


, 


|IM(X, Y)? 2 N(1-oaYP +cX +4X? +..), (4) 


the so-called linear representation, where a is a complex 
number. The real part of o gives the linear term in Y 
for the Dalitz plot density, a = 2#(a), and the quadratic 
term is b = R(a)? + S(a)?, where R(a) and S(a) are 
the real and imaginary parts of a, respectively. The two 
representations are equivalent if b > a?/4, i.e., b should 
be at least larger than zero. Therefore, a negative value 
for b demonstrates that the ansatz of Eq. (4) does not 
describe the data. 

Experimentally, the decays of the y! — nrm and 
n! — nn°n® have only been explored with limited statis- 
tics so far. The matrix elements for n — nnr have 
been studied by the CLEO (using only Eq. (4)) [8], 
VES [9] and BESIII [10] Collaborations. The most re- 
cent measurement of n! + nnr? is from the GAMS-47 
experiment [11], complementing older results reported by 
the GAMS-2000 Collaboration [12]. Discrepancies in the 
Dalitz plot parameters both for the charged and neutral 
decay channels are obvious from those experiments. 

In addition, the Dalitz plot for 7’ — nr°r° is expected 
to be affected by a cusp due to the ^^ mass threshold. 
The size of this effect is predicted to be about 8% [7]. A 
analogous cusp has been observed in K * — «^ «9s? [13] 
decay and allows to determinate the v7 S-wave scattering 
lengths [14-17]. 

The dynamics of the decays 7! — nata7~ and n > 
r9? is studied in this work using n’ mesons produced in 
J/ — yn! decay. The present data sample of 1.31 x 10? 
J/w events accumulated with the BESIII detector is 
about 5 times of that used in the previous BESIII anal- 
ysis [10]. 


II. DETECTOR AND MONTE CARLO 
SIMULATION 


The BESIII detector is a general-purpose magnetic 
spectrometer with a geometrical acceptance of 93% of 
47, and is described in detail in Ref. [18]. The detec- 
tor is composed of a helium-based drift chamber (MDC), 
a plastic-scintillator time-of-flight system (TOF), and a 
CsI(T1) electromagnetic calorimeter (EMC), all enclosed 
in a superconducting solenoidal magnet providing a 1.0 T 
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magnetic field (0.9 T in 2012). The solenoid is sup- 
ported by an octagonal flux-return yoke with resistive- 
plate counters interleaved with steel for muon identifica- 
tion (MUC). 

The GEANTA [19] based Monte Carlo (MC) simula- 
tion software package BOOST [20] describes geometry 
and material of the BESIII detector, as well as detector 
response. It is used to optimize the event selection cri- 
teria, estimate backgrounds and determine the detection 
efficiencies. The production of the J/4» resonance is sim- 
ulated with KKMC [21, 22], while the decays are gener- 
ated with EVTGEN [23, 24] for established modes using 
world-average branching fractions [1], and by LUND- 
CHARM [25] for the remaining decays. An inclusive 
MC sample of 1.2 x 10° J/z events is used to study the 
potential background contributions. The analysis is per- 
formed in the framework of the BESIII offline software 
system (BOSS) [26]. 


II. MEASUREMENT OF THE MATRIX 
ELEMENTS FOR 7 > nati 


For the reconstruction of J/4» — yn! with n! — neta 
and 7 — yy, candidate events must contain two tracks 
with opposite charge and at least three photons. Each 
charged track reconstructed from the MDC hits is re- 
quired to have a polar angle in the range | cos0| < 0.93 
and to pass the interaction point within +10 cm along 
the beam direction and within +1 cm in the plane per- 
pendicular to the beam. Photon candidates are recon- 
structed using isolated clusters of energy deposited in 
the EMC and required to have a deposited energy larger 
than 25 MeV in the barrel region (|cos0| < 0.80) or 
50 MeV in the endcap region (0.86 < |cos0| < 0.92). 
The energy deposited in nearby TOF counters is included 
to improve the reconstruction efficiency and energy res- 
olution. To eliminate clusters associated with charged 
tracks, the angle between the photon candidate and the 
extrapolation of any charged track to the EMC must be 
larger than 10?. A requirement on EMC cluster timing 
with respect to the event start time (0 < T < 700 ns) 
is used to suppress electronic noise and energy deposits 
unrelated to the event. 

Since the radiative photon from J/w decay is always 
more energetic than those from 7 decay, the photon can- 
didate with the maximum energy in an event is assumed 
to be the radiative one. For each c *- ^yy combina- 
tion, a six-constraint (6C) kinematic fit is applied and 
the X2c is required to be less than 100. The fit enforces 
energy-momentum conservation and constrains the in- 
variant masses of yy and rx ^. to the nominal 7 and y’ 
masses, respectively. If there are more than three selected 
photons in an event, the combination with the smallest 
Xéc is retained. 

To estimate the background contribution, an alterna- 
tive data sample is selected without applying the 7 and 
n’ mass constraints in the kinematic fit. The zz yy in- 


Entries/(2.5MeV/c?) 


085 09 095 1 105 11 
M(x*rn) (GeV/c?) 


Figure 1. Invariant mass spectrum of 7^ 7 candidates with- 
out 7 and 7’ mass constraints applied in the kinematic fit, and 
requiring the yy invariant mass within the 7 signal region. 


variant mass spectrum is shown in Fig. 1 after requiring 
the yy invariant mass within the 77 signal region, (0.518, 
0.578) GeV/c?. A clear 1 signal is observed with a low 
background level. In addition, a sample of 1.2 x 10? 
inclusive MC J/w decays is used to investigate poten- 
tial backgrounds. Using the same selection criteria for 
the MC sample, no peaking background remains around 
the 7 signal region. From this MC sample, the back- 
ground contamination is estimated to be about 0.396. 
This is consistent with an estimation obtained from an 
unbinned maximum likelihood fit to the M (rm) dis- 
tribution, where the signal is described by the MC sim- 
ulated shape convoluted with a Gaussian function rep- 
resenting the resolution difference between data and MC 
simulation, and the background contribution is described 
by a third-order polynomial function. We therefore ne- 
glect the background contribution in the determination 
of the Dalitz plot parameters. 


After the above requirements, 351,016 1 — nrm” 
candidate events are selected, with a background contri- 
bution of less than 0.396. Figure 2 shows the Dalitz plot 
in the variables X and Y for the selected events. The 
corresponding projections on .X and Y are shown as the 
dots with error bars in Figs. 3 (a) and (b), respectively. 


Unbinned maximum likelihood fits to the data are per- 
formed to determine the free parameters in the decay am- 
plitude squared (Eqs. (3) and (4)). To account for the 
resolution and detection efficiency, the amplitude squared 
is convoluted with a function o(X, Y) parameterizing the 
resolution, and multiplied by a function e(.X, Y) param- 
eterizing the detection efficiency. Both functions are de- 
rived from MC simulations. Two double Gaussian func- 
tions are used for c(X,Y), while £(X,Y) is estimated 
as the average efficiencies of local bins. With the nor- 
malization, one derives the probability density function 
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Figure 2. Dalitz plot for n > nat from data. 


P(X,Y), which is applied in the fit: 


POLY) = 
—IM(XY)P@o(KY) XY) 5) 
Top (MX, YP @ o(X, Y) -e(X, Y))dXdY 


The integral over the full Dalitz plot range (DP) gives 
the normalization factor in the denominator. The fit is 
done by minimizing the negative log-likelihood value 


Nevent 
-n£=- M mP(Xi,%), (6) 


i=l 


where P(X;, Yi) is evaluated for event i, and the sum 
runs over all accepted events. 

Imposing charge conjugation invariance by setting the 
coefficient of odd powers in X (c) to zero in the general 
representation, the fit yields following parameters: 


a = —0.056 + 0.004, 
b = —0.049 + 0.006, (7) 
d = —0.063 + 0.004. 


Here, the uncertainties are statistical only. The corre- 
sponding correlation matrix of the fit parameters is 


a | —0.417 —0.239 | . (8) 
b 0.292 


Projections of the fit result on X and Y are illustrated 
as the solid histograms in Fig. 3. 

To check for the existence of a charge conjugation vio- 
lating term, an alternative fit with parameter c left free 
is performed. The resultant value, c = (2.7+2.4) x 1073, 
is consistent with zero. Compared with the nominal fit 
results, the parameters a, b and d are almost unchanged, 
and the statistical significance for a non-zero value of pa- 
rameter c is determined to be 0.76 only. 


Alternative fits including the extra terms fY? --gX?Y 
or eXY + hXY? - 1X? in the general representation are 
also performed, resulting in f = —0.004 + 0.012, g = 
0.008 + 0.010 or e = 0.005 +0.007, h = 0.004 0.006, and 
l = 0.007+0.013, respectively, while the other parameters 
are unchanged. 

A fit based on the linear representation is also per- 
formed and yields the following values: 


(a) = —0.034 + 0.002, 
S(a) = 0.000 + 0.019, (9) 
d = —0.053+4 0.004. 


The imaginary part of a is consistent with zero. This can 
be understood by the observation that the coefficient b 
in the general representation is negative. 

Subsequently we will consider the fit result with S(a) 
fixed at zero. The parameters * (o) and d and their un- 
certainties remain the same as in Eq. (9), and the cor- 
relation coefficient between R(a) and d is —0.137. The 
log-likelihood value is lower by 33.9 compared with the fit 
using the general representation, which indicates the in- 
compatibility of the linear representation with the data. 
Projections on .X and Y based on this result are illus- 
trated as the dashed histograms in Fig. 3. The presented 
residuals show that the fit is slightly worse to describe 
the data in Y projection comparing to the general one. 

The potential charge conjugation violating is also 
checked in the linear representation by performing an 
alternative fit with a free parameter c. The resultant 
value, c = (2.7+2.4) x 107, is also consistent with zero, 
while the parameters (o) and d are almost unchanged 
compared with the nominal fit results. 


IV. MEASUREMENT OF THE MATRIX 
ELEMENT FOR THE DECAY ;/ > n7°r® 


In the reconstruction of J/4» — yn! with n! — nz? 
and n/nx° — yy, candidate events must have at least 
seven photons and no charged track. The selection cri- 
teria for photon candidates are the same as those for 
n! — nnr, except that the requirement on the angle 
between photon candidates and any charged track is not 
used. A requirement of EMC cluster timing with respect 
to the most energetic photon (—500 < T < 500 ns) is 
also used. The photon with the largest energy in the 
event is assumed to be the radiative photon originating 
from the J/w decay. For the remaining clusters, pairs of 
photons are combined into 2/1 — yy candidates, which 
are subjected to a one-constraint (1C) kinematic fit by 
constraining the invariant mass of the photon pair to be 
the nominal 7° or 7 mass. The x? for this 1C kinematic 
fit is required to be less than 25. To suppress 7? miscom- 
binations, the 7° decay angle Üaccay, defined as the polar 
angle of one of the decay photons in the yy rest frame 
with respect to the 7° flight direction, is required to sat- 
isfy | cos Ogecay| < 0.95. Then an eight-constraint (8C) 
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Figure 3. Projections of the fit results onto (a) X and (b) Y in the general (solid histograms) and linear (dashed histograms) 
representations for n — rz *-- , where the dots with error bars represent data. 
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Figure 4. Invariant mass spectrum of 7°7r°n candidates with- 


out 7 and 7’ mass constraints applied in the kinematic fit, and 
requiring the yy invariant mass within the 7 signal region. 


kinematic fit is performed for the ynr°r° combination 


enforcing energy-momentum conservation and constrain- 
ing the invariant masses of the three photon pairs and the 
nr’? combination to the nominal 7°/7 and r masses. 
If more than one combination is found in an event, only 
the one with the smallest x2. is retained. Events with 
Xc < 100 are accepted for further analysis. 

To estimate the backgrounds, an alternative selection 
is performed where 7 and m' mass constraints in the 
kinematic fit are removed. The resulting «95975 invari- 
ant mass spectrum is shown in Fig. 4, after requiring 
the yy invariant mass within the 7 signal region, (0.518, 


0.578) GeV/c?. The inclusive MC study shows that 
the surviving backgrounds mainly consist of the peaking 
background 7 — 7077? and a flat contribution from 
J/A — wn with w > yr? and n — n?n9«?. From this 
MC sample, the background contamination is estimated 
to be about 0.996, which is consistent with the estima- 
tion obtained from a fit to M(nr°x°), and therefore ne- 
glected in the determination of the Dalitz plot parame- 
ters. In the fit, the signal is described by the MC sim- 
ulated shape convoluted with a Gaussian function rep- 
resenting the difference of mass resolution between data 
and MC simulation. The shape and the yield of peak- 
ing background n — «9595? are fixed according to the 
dedicated MC simulation [27]. A third-order polynomial 
function is used to represent the smooth background con- 
tribution. 

After the above requirements, 56,249 n > nr°r° can- 
didate events are selected, with a 0.9% background level. 
The Dalitz plot of selected events is displayed in Fig. 5. 
The corresponding projections on X and Y are shown as 
the dots with error bars in Fig. 6. 

As in the analysis of the n” — nrm, an unbinned 
maximum likelihood fit method is used to determine the 
matrix elements. The resolution is described with two 
double Gaussian functions, and the detection efficiencies 
in different X and Y bins are obtained from MC simu- 
lation. From a dedicated study with the control sample 
of J/ip — nrn’, we find that the reconstruction effi- 
ciency for 7° candidate differs significantly between data 
and the MC simulation at low 7° momenta. Thus, to 
describe the detection efficiency more accurately, an ef- 
ficiency correction depending on the 7° momentum is 
carried out, and the error of this correction will be con- 
sidered in systematic uncertainty. 
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Figure 5. Dalitz plot for 7’ — nx°x° from data. 


Considering the strict constraint from the symmetry 
of the wave function, only the fits without odd powers of 
X are performed. The fit based on the general represen- 
tation yields the coefficient (with statistical uncertainties 
only) and the corresponding correlation matrix: 


a = —0.087 + 0.009, 
b = —0.073 + 0.014, (10) 
d = —0.074 + 0.009, 


a | 0.495 —0.273 | . (11) 
b 0.273 


Similarly to the case 7’ — natn, the fit gives a negative 
value of the coefficient b. The projections of the fit results 
on X and Y are shown as the solid histograms in Fig. 6. 

Extra terms fY? and gX?Y in the general represen- 
tation are also added in an alternative fit, resulting in 
f = —0.023 + 0.028 and g = 0.024 + 0.025. The signifi- 
cance for non-zero values of f and g is 0.66. 

In the fit based on the linear representation, the imag- 
inary part of a (fitted to be 0.000 + 0.038) does not con- 
tribute to the fit quality, as in the y' — nat case. 
Thus, the nominal fit omitting S(a) gives the results 
(statistical uncertainties only): 


(a) = —0.054 + 0.004, 
d — 0.061 + 0.009. 


(12) 


'The correlation coefficient between the two parameters 
is —0.170. Compared to the fit based on the general rep- 
resentation, the log-likelihood value is reduced by 13.7. 
Projections on X and Y are illustrated as the dashed 
histograms in Fig. 6. Again, the fits based on the two 
different representations give similar results for the X 
projections, but slightly worse for Y in the linear case. 


V. SYSTEMATIC UNCERTAINTIES 


Various sources of systematic uncertainties on the mea- 
sured Dalitz plot matrix elements have been investigated, 
including tracking efficiency, kinematic fit, efficiency cor- 
rection, and resolution. For the decay n’ — nn??, addi- 
tional uncertainties associated with photon miscombina- 
tion, 7° and 7 reconstruction are also considered. 

Differences between data and MC simulation for the 
tracking efficiency of charged pions are investigated using 
the control sample J/q > ppr*- . A momentum de- 
pendent correction on the detection efficiency is obtained 
by comparing the efficiency between data and MC simu- 
lation. Similarly, a momentum dependent correction for 
the 7 reconstruction efficiency is obtained with the con- 
trol sample of J/y — -ygm*m«- . Then alternative fits 
are performed by incorporating the efficiency corrections 
for charged pions or y. Changes of the matrix elements 
with respect to the nominal results are assigned as the 
systematic uncertainties. A momentum-dependent 7° re- 
construction efficiency correction has been applied in the 
nominal fit; the associated systematic uncertainties are 
estimated by changing the correction factor by one of its 
standard deviation and repeating the fit. 

'The possible miscombination of photons in signal MC 
samples has been studied by matching the generated pho- 
ton pairs to the selected 7° or 7 candidates. The frac- 
tion of events with wrong combinations is determined 
to be 2.796 for n' — nr°n®. Alternative fits are per- 
formed to the MC simulated sample with only truth- 
tagged events and the ones including miscombinations, 
individually. The difference between those two results 
are taken as the the systematic uncertainties. 

'To estimate the uncertainties associated with the kine- 
matic fitting procedure, the fit results are compared us- 
ing a 4C (6C) instead of 6C (8C) kinematic fit for n' > 
qrmn- (m + nr°r°) and the corresponding changes in 
the fit parameters are taken as systematic uncertainties. 

To estimate the uncertainties associated with the ef- 
ficiency correction in Eq. 5, we change the Dalitz plot 
variables X and Y to the so-called square Dalitz plot vari- 
ables M (rr)? and cos 0, where 0 is the angle between the 
two pions in the rest frame of nm. Alternative fits are per- 
formed with the efficiency correction based on the newly 
defined Dalitz plot variable and the resultant changes of 
the matrix elements with respect to the nominal results 
are assigned as systematic uncertainties. 

To estimate the impact from the non-flat resolution in 
the X —Y plane, the biases from input/output checks are 
taken as the systematic uncertainties. The impact from 
different resolutions of the Dalitz plot variables between 
data and MC simulation is estimated by alternative fits 
varying the resolutions by +10%. It is found that the 
change of the results is negligible. The effect of neglect- 
ing the residual background is checked by alternative fits 
including MC simulated backgrounds and found to be 
insignificant. 

All of the above uncertainties are summarized in Ta- 
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Figure 6. Projections of the fit results on (a) X and (b) Y in general (solid histograms) and linear (dashed histograms) 
representations for 7’ — qz ^x?, where the dots with error bars represent data. 


ble I. Assuming all the sources of systematic uncertainty 
are independent, the total systematic uncertainties for 
the matrix elements are obtained by adding the individ- 
ual values in quadrature, shown in the last row of Table I. 


VI. SEARCH FOR CUSP EFFECT IN y > zs? 


After the event selection presented in section IV, the 
distributions of the square of the 7°7° invariant mass, 
M(n°n°)?, for data and phase space distributed MC sim- 
ulation are depicted in Fig. 7 (a), together with the Dalitz 
plot fit result based on the general representation. To 
check the cusp effect, Fig. 7 (b) shows the M (1959)? dis- 
tribution of data divided by the phase space distribution. 
We observe no statistically significant structure near the 
ntn mass threshold, which would indicate a cusp. The 
prediction in Ref. [7] was based on the matrix elements 
measured by VES experiment [9]. 


VII. SUMMARY 


With a sample of 1.31 x 10? J/w events collected 
with the BESIII detector, clean samples of 351,016 7’ > 
nit events and 56,249 y — nz? events are selected 
from J/14b radiative decays. Then the most precise mea- 
surements of the matrix elements for the n! > qgr*«- 
and 7’ — nnz?«? decays as well as a search for the cusp 
effect in n! + n7°7r° are performed. 

Both the general and the linear representations are 
used to determine the matrix elements and the corre- 
sponding results are summarized in Table II including 


the systematic uncertainties. The matrix elements for 
both decays are in reasonable agreement and more pre- 
cise than the previous measurements [8-11]. The results 
for n! + nnr supersede the previous BESIII measure- 
ment [10], which used a subsample of the present data. 
As reported in Ref. [10], the discrepancy of the parame- 
ter a for n! + nx*+n~ with respect to the VES value [9] is 
evident, which, at present, stands at about 3.8 standard 
deviations. The values of parameter c in n — nat are 
all consistent with zero within one standard deviation in 
both representations, in agreement with charge conjuga- 
tion conservation in the strong interaction. In addition, 
a discrepancy of 2.6 standard deviations for parameter a 
is observed between 7! — nrm and 7! — nnr? pro- 
cesses, indicating isospin violation. However, the result is 
not statistical significant enough to firmly establish such 
a violation, and additional effects, e.g., radiative correc- 
tions [7], should be considered in future experimental and 
theoretical studies. 

A comparison between the results obtained from gen- 
eral representation and the theoretical predictions within 
the framework of U(3) chiral effective field theory (EFT) 
incorporating with a relativistic coupled-channels ap- 
proach [4] is given in Table IL. In general, our results 
are compatible with the theoretical expectations. How- 
ever, the theoretical prediction for parameter a from 
n! —^ natn is about two times larger than our result and 
the discrepancies on parameter d for both n > qymr*«- 
and n! — nr°r® are about 4 standard deviations. Re- 
cent calculations with large-Nc ChPT at next-to-next- 
to-leading order [5] are compatible with our measurement 
of parameter a due to the large theoretical uncertainty. 

As previously mentioned, the linear and general repre- 
sentations are equivalent for the case of b > a?/4. How- 


Table I. Relative systematic uncertainties of the matrix elements in the generalized and linear representations. 
n — nnn 


Source General Rep. Linear Rep. General Rep. Linear Rep. 
a b d | (o) d a b d R(a) d 


Tracking efficiency — |0.0018 0.0044 0.0021/0.0015 0.0013 
n? efficiency 


0.0006 0.0007 0.0001|0.0004 0.0003 
0.0012 0.0014 0.0001}0.0005 0.0003 
sis TA re Pe .. [0.0002 0.0024 0.0013 ]|0.0004 0.0009 
0.0009 0.0035 0.0024|0.0007 0.0031 0.0041 0.0031 0.0020}0.0005 0.0016 
Efficiency presentation 10.0009 0.0002 0.0007 0.0005 0.0007 0.0002 0.0005 0.0004 :0.0002 0.0051 


n efficiency 
Photon miscombination 


Kinematic fit 


Resolution 0.0006 0.0093 0.0004|0.0005 0.0015 |0.0044 0.0021 0.0030}0.0005 0.0049 
Total 0.0023 0.0057 0.0034|0.0018 0.0038 }0.0061 0.0047 0.0038|0.0010 0.0052 
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Figure 7. (a) The distribution of the n?n’ invariant mass squared. In the upper plot, the dots with error bars represent data, 
the dashed histogram is from phase space distributed MC sample, and the solid histogram is the projection of the Dalitz plot fit 
in general representation; the lower gives the distribution of the ratio between data and the fit projection. (b) The M(1°n°)? 
distribution of data divided by phase space distribution. The vertical line is the v * c mass threshold. 


Table II. Experimental and theoretical values of the matrix element parameters for n! — gx *«- and n > nr. The values 
for parameter c and S(q) are given for comparison with previous experiments. 


Para- y = ain 
meter Large No [5] VES [9] This work EFT [4] GAMS-4z [11] This work 
a E 0.011 —0.098 4 t 0.018 —0.056 4 0.127 + 0.009 —0.067 + 0.016 —0.087 + 0.009 + 0.006 
b + 0.034 —0.050 4 - 0.032 —0.049 4 0.049 + 0.036 —0.064 + 0.029 —0.073 + 0.014 + 0.005 
c E - 0.018 (2.7 + 2.4 + 1.8) x 10? a m "s 
d E 0.019 —0.092 4 -0.019 —0.063 4 E 0.004 |-+0.011 + 0.021 —0.067 + 0.020 —0.074 + 0.009 + 0.004 
R(a) » T —0.072 + 0.014 —0.034 + 0.002 + 0.002 e —0.042 + 0.008 —0.054 + 0.004 + 0.001 
S(o) T T 0.000 2- 0.100 — 0.000 + 0.019 + 0.001 = 0.000 + 0.070 0.000 + 0.038 + 0.002 
c +0.020 + 0.019 (2.7 4: 2.4 + 1.5) x 10? a E EN 
d —0.066 + 0.034 —0.053 + 0.004 + 0.004 Bs —0.054 + 0.019 —0.061 + 0.009 + 0.006 
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ever, the coefficients b for the Y? term are negative with 
5.8 and 4.9 standard deviations to zero for 7! — nnr" 
and 7! — nxn, respectively, which implies that these 
two representations can not provide an identical descrip- 
tion of data. In case of the linear representation, the re- 
sults are in agreement with previous measurements and 
also provide a reasonable description on .X projection 
for both decay modes. However, the goodness of fit on 
the Y projections are worse than the general one. This 
is consistent with the conclusion reported by the VES 
collaboration [9] that the linear representation can not 
describe the data well. 

We also attempt to search for the cusp effect in the 
decay n — nr?7?. Inspection of the 195? mass spectrum 
around the 7*7 mass threshold does not show evidence 
of a cusp with current statistics. 
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